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A Toxoplasma dense granule protein, GRA24, 

modulates the early immune response 

to infection by promoting a direct 

and sustained host p38 MAPK activation 
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Toxoplasma gondii, the causative agent of toxoplasmosis, is an obligate intracellular proto- 

p ^ zoan parasite that resides inside a parasitophorous vacuole. During infection, Toxoplasma 

X actively remodels the transcriptome of its hosting cells with profound and coupled impact 

W on the host immune response. We report that Toxoplasma secretes GRA24, a novel dense 

^ granule protein which traffics from the vacuole to the host cell nucleus. Once released into 

I the host cell, GRA24 has the unique ability to trigger prolonged autophosphorylation and 

^ nuclear translocation of the host cell p38a MAP kinase. This noncanonical kinetics of p38a 

U activation correlates with the up-regulation of the transcription factors Egr-1 and c-Fos 

Q and the correlated synthesis of key proinflammatory cytokines, including interleukin-1 2 

H-j and the chemokine MCP-1, both known to control early parasite replication in vivo. Re- 

^ markably, the GRA24-p38a complex is defined by peculiar structural features and uncovers 

I a new regulatory signaling path distinct from the MAPK signaling cascade and otherwise 

commonly activated by stress-related stimuli or various intracellular microbes. 
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Toxoplasma gondii is a widespread obligate intra- 
cellular protozoan parasite causing toxoplasmosis, 
a potentially severe disease in imniunocompro- 
niised or congenitaUy infected humans. Once 
it has invaded its hosting cell, the tachyzoite 
grows inside a parasitophorous vacuole (PV) 
in any type of nucleated cell, where it directs 
profound changes in the host cell transcriptome 
(Martin et al., 2007; Plattner and Soldati-Favre, 
2008), which is known to tightly precede the 
host cell proteome (Nelson et al., 2008). In this 
regard, when the hosting cells are of mouse 
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origin, IL-12 and IFN-7 were recognized as 
a pivotal immune axis contributing not only 
to the control of the tachyzoite population 
size (Sher et al., 2003; Mashayekhi et al, 2011; 
Goldszmid et al., 2012) but also to the T. gondii 
tachyzoite-bradyzoite developmental transition. 
Interestingly, once delivered within the host 
cell cytoplasm, lineage-specific polymorphic 
tachyzoite molecules have been shown to coun- 
teract this IL-12p70— IFN-7 axis, either acting 
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□ Signal peptide 

■ Nuclear Localization Signal 



■ Internal repeat 
□ KIM/D-motif 



□ S/TP target site 
* STOPcodon 



ME49_030180 MLQMARYTANICAVSICSLVLWALSVGILPTPDWKDRVKMGGTESGPFV 
VEG_0 25320 MLQMARYTANICAVSICSLVLWALSVGILPTPDWKDRVKMGGTESGPFV 
GT1_117710 MLQMARYTVNICAVSICSLVLWALSVDILPTPDWKDRMKMGGTESGPFV 

ME49_030180 LQVCASDPLLHDPKERESGSDSTRGYHGGSSSGGSSSRQGTTVRSDAGLS 
VEG_025320 LQVCASDPLLHDPKERESGSDSTRGYHGGSSSGGSSSRQGTTVRSDAGLS 
GT1_117710 LQVCASDPLLHAPKERESGSDSTRGYHGGSSSGGSSSRQGTTVRSDAGPS 

ME4 9_03 0 180 SQSSQSSASTSAKTSENHRQGPAFLTSVFRKGETPALDWVPYGTLEGAKW 
VEG_025320 SQSSQS S ASTS AKTSENHRQGPAFLTS VFRKGETPALDWVPYGTLEGAKW 
GT1_117710 SQSSQS S ASTS AKTSEKHQQGPAFLTS VFRKGETPALHWVPYGTLEGAKW 

ME49_030180 HPGQQKSKRRSSATTSRRQGASHSGNPGQLPAPRGGLQPTTTLSGTGGQP 
VEG_025320 HPGQQKSKRRSSATTSRRQGASHSGNPGQLPAPRGGLQPTTTLSGTGGQP 
GT1_117710 HPGQQKSKRRSSATTSRQQGASHSGNPGQLPAPRGGLQPTTTLSGTAGQP 

ME4 9_03 0 180 RTDSTDEGAAATSVIPNRSGYPQPVPYLIHPVGFLSGDYRSLGMLGLVPS 
VEG_025320 RTDSTDEGAAATSVIPNRSGYPQPVPYLIHPVGFLSGDYRSLGMLGLVPS 
GT1_117710 RTDSTDEGAAATSVIPNRSGDPQPVPYLIHPVGFLSGDYNSLGHSGLVPS 
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Figure 1. The gene GRA24 generates splice 
variants encoding two distinct proteins. (A) Sche- 
matic representation of the 6R/A24 alternative splice 
variants and their corresponding GRA24 protein iso- 
forms. The putative signal peptide is shown in yellow 
(the cleavage site is predicted by SignalP). A putative 
nuclear localization signal is in green (as predicted by 
PSORT). The KIM/D-motif (in red square) and the 
internal repeats (colored in red) are shown. (B) GRA24 
amino acid alignment. The predicted amino acid se- 
quence for the primary translation product of the 
GRA24 gene is shown for type I (TGGJ1_ 1 17710], 
II {TGME49_030180), and III (TGVEG_025320) para- 
sites. Single amino acid polymorphisms are indicated 
for type II and III (in blue) versus type I (in red). The 
alignment was done by ClustalW2 (EMBL-EBI). 
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directly — e.g., disruption of the autonomous host cell protecting 
immune processes (Howard et al., 2011) — or indirectly — e.g., 
subversion of the transcriptonie machinery (Hunter and Sibley, 
2012). The majority of these tachyzoite molecules identi- 
fied so far are transiently released early on during host cell 
invasion from the rhoptry secretory organelles (Hunter and 
Sibley, 2012). Chief among them, the polymorphic rhoptry 
protein ROP16 directs tyrosine phosphorylation of STAT3 
and STAT6 (Saeij et al, 2007;Yamamoto et al, 2009; Butcher 
et al., 2011). Once enclosed in the PV, additional effectors, 
including GRA15 (Rosowski et al., 2011) and GRA16 
(Bougdour et al., 2013), are secreted from dense granules 
(DGs) and traffic to the host cell. Notably, GRA16 once ex- 
ported beyond the PV membrane (PVM) reaches the host 
cell nucleus, where it positively modulates genes involved in 
cell cycle progression and the p53 tumor suppressor pathway 
(Bougdour etal.,2013). 



Considerable effort is currently being made to unravel 
the intricacy of the parasite and host interplay that leads, for 
the medically relevant type II and III strains, to the control of 
acute parasite multiplication and dissemination yet allows 
persistence of cryptic parasites in deep tissues. What is known 
is the need for a strong and persistent Thl response and 
the correlated and marked proinflammatory cytokine IL-12 
production by innate immune cells, including macrophages 
(M0), CD8a^ dendritic cells, and neutrophils (Sher et al., 
2003; Mashayekhi et al, 2011; Goldszmid et al., 2012). IL-12, 
in turn, activates natural killer and T cells to produce IFN-7, 
which triggers parasite killing through the activation of IRGs 
(IFN-regulated GTPases; Howard et al., 2011). In mice lacking 
IL-12 or IFN-"y, parasite expansion throughout the host is 
not controlled, causing host death (Gazzinelli et al., 1994; 
Scharton-Kersten et al., 1996;Yap et al, 2006). Although IFN-7 
is the main mediator of resistance to Toxoplasma infection. 
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IL-12 synthesis appears as the key cytokine to be tightly regu- 
lated in a strain-specific manner by Toxoplasma. In this regard, 
macrophages are classically (or Ml) activated by a type II 
strain while they are alternatively (or M2) activated by hyper- 
virulent type I and avirulent type III strains Qensen et al., 
2011).While GPJ\.15 type II allele contributes to Ml activation 
by activating NF-kB and downstream proinflammatory cyto- 
kines, in particular IL-12p70 (Jensen et al., 2011; Rosowski 
et al, 2011), ROP16 triggers a STAT3/6-dependent M2 re- 
sponse that prevents IL-12p70 from being synthesized (Saeij 
et al, 2007; Butcher et al, 2011;Jensen et al, 2011). 

Sustained IL-12p70 production in tissues hosting T. gondii is 
not exclusively controlled by NF-kB but could also require the 
MAPK pathways (Kim et al, 2005; Masek et al, 2006). MAPK 
pathways encompass a series of kinases that activate each other 
through an orchestrated cascade (MAP4K, then MAP3K, then 
MAP2K) that ends with the phosphorylation of a specific 
MAPK that positively or negatively regulates the expression of 
suites of genes by activating effector proteins, in particular tran- 
scription factors. A wide range of stimuli activates MAPKs, in- 
cluding inflammatory cytokines, osmotic stress, LPS, growth 
factors, osmotic shock, and UV light (Krishna and Naraiig, 
2008; Cuadrado and Nebreda, 20 10). Although it has been pro- 
posed that Toxoplasma bypasses the MKK-dependent MAPK 
phosphorylation cascade to induce autophosphorylation of the 
p38ct MAPK and to promote IL-12p70 production (Kim et al, 
2005), the tachyzoite molecule subverting the host p38a 
MAPK pathway remains unknown. 

Within this context, we characterize a novel DG-like pro- 
tein, hereafter referred to as GRA24, that is exported beyond 
the PV to the host cell nucleus. We discovered that GRA24 
directly interacts with p3Ba, leading to an unusually persis- 
tent autophosphorylation and activation of the host kinase. 
3D modeling of the complex predicts unusual structural fea- 
tures of the complex that highlight a novel pathway for p38a 
extended activation. The GRA24-dependent p38a activation 
was documented in both human stromal cells and mouse 
macrophages, the latter up-regulating the expression of mul- 
tiple proinflammatory cytokines (i.e., IL-12p40 and MCP-1/ 
CCL2) in a strain-specific manner. Our results demonstrate a 
substantial role for GRA24 in the control of parasite burden 
in the gut of infected mice. Finally, as a link between the cyto- 
kine response and GRA24-p38a activation, we identified sev- 
eral transcription factors up-regulated in a GRA24-dependent 
fashion, among which are Egr-1 and c-Fos. Collectively, our 
findings identify a GRA24-dependent sustained p38a kinase 
activation pathway that could contribute shaping and modulat- 
ing immune response — at least the one driven by the T. gondii 
tachyzoite developmental stage. 

RESULTS 

GRA24 is a novel DG-like protein exported 
to the host cell nucleus 

The gene TGME49_030180 encoding GRA24 was originally 
identified in srlico together with GRA16 (Bougdour et al., 
2013) as strong parasite candidate genes with attributes to 
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Figure 2. GRA24 is a novel DG protein exported beyond the PV 
into the host cell nucleus. (A) Time course of GRA24 secretion and 
export to tine Inost cell nucleus. Human fibroblasts (HFF) were infected 
with parasites expressing an HAFlag (HF)-tagged copy of GRA24 
{RHI(u80GRA24-Hf], fixed, and stained with anti-HA antibodies (red). 
Images are representative of at least four experiments. Bar, 20 \im. 
(B) GRA24-HAFIag is contained in cytoplasmic organelles distinct from the 
apical micronemes and rhoptries and partially colocalizing with DG pro- 
teins. (B-E) Free parasites were fixed, permeabilized, and successively 
stained with antibodies raised against HA to label GRA24 (green) and then 
with antibodies (red) against the micronemal protein MIC2 (B), the rhoptry 
protein toxofilin (C), and the DG proteins GRA7 (D) and GRAl (E). The top 
panels show the DIC images corresponding to the IFA (bottom panels), and 
the right frames show 3D images processed using MetaMorph and Imaris 
software. Images are representative of at least two experiments. Bars, 5 pm. 
(F) IFA of GRA24-HF expressed by the strong promoter of GRAl. HFF 
monolayers were infected with RH<ruSO parasites expressing ectopically 
pGRAl-GRA24-HF 18 h after infection, cells were fixed and secretion of 
GRA24 was monitored using anti-HA antibodies (green). Images are repre- 
sentative of at least two experiments. Bar, 10 pm. (G) Higher magnification 
of a vacuole showing GRA24-HF signal in the vacuolar space. Bar, 2.5 pm. 
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Figure 3. GRA24 forms a dimeric com- 
plex with the p38o( MAPK. (A) Immuno- 
fluorescence analysis (IFA) of GRA24-HAFIag 
ectopically and stably expressed in 293-TRex 
cell line. Cells were either left untreated (— ) 
or treated with 1 pg/ml tetracycline for 20 h 
before fixation and staining with anti-HA 
antibodies (red) and Hoechst DNA-specific dye 
(blue). Bar, 2.5 pm. (B) GRA24-associated 
polypeptides were purified from nuclear ex- 
tracts of 293-TRex cells tetracycline-induced 
to express GRA24-HF. Size exclusion chroma- 
tography (SEC) of GRA24-containing com- 
plexes after Flag affinity selection. Fractions 
were analyzed by silver staining and immuno- 
blotting to detect GRA24-HF (anti-HA), total 
p38c«, and Thr180/Tyr182 phosphorylated 
*P-p38a. Data are representative of two ex- 
periments. (C) Mass spectrometry analysis of 
SEC fraction 22. Identity of the proteins with 
their respective number of peptides 
is indicated on the right. (D) GRA24-p38oi 
association detected in HFF and J774 M0 
infected by RHku80 GRA24-HF (18 h). Tagged 
proteins were Flag-immunoprecipitated from 
extracts and eluates were analyzed by immuno- 
blotting. Data are representative of two ex- 
periments. (E) GRA24-associated polypeptides 
were purified from nuclear extracts of 293- 
TRex cells induced to express an N-terminal 
HAFIag-tagged GRA24. SEC Fractions were 
analyzed by immunoblotting to detect 
HAFIag-6RA24 (anti-HA), total p38a, and 
Thr180/Tyr182 phosphorylated *P-p38a. Data 
are representative of two experiments. 



both target the parasite secretory pathway and to reach the 
nuclei of T. gondii tachyzoite-hosting cells. GRA24 encodes 
a 542-aa protein with a predicted N-terniinal signal peptide, 
an internal putative bipartite nuclear localization sequence, 
and two repeats, Rl and R2, at the C-terminal end (Fig. 1 A). 
GRA24 protein is unique because it has no significant ho- 
mology with any proteins, not even with the close relative 
Neospom caninum proteins. In Toxoplasma, GRA24 harbors 
significant polymorphism between the allele shared by types 
II and III clonal lineages and the allele typifying type I clonal 
lineage (Fig. 1 B). In fibroblasts hosting tachyzoites express- 
ing a carboxy-terminal HAFIag-tagged version of GmV24, 
the protein is distinctly detected in the host cell nucleoplasm 
at 16 h after invasion and later, whereas there is no detectable 
protein left in the intracellular tachyzoites (Fig. 2 A). In con- 
trast, in free tachyzoites, GRA24-HAFlag delineates punc- 
tuate structures throughout the parasite cytoplasm which 
occasionally overlap with the canonical DG proteins GRAl 
and GRA7 (Fig. 2, D and E) while being excluded from the 
two typical apical secretory organelles using toxofilin and 
MIC2 as markers for rhoptry and niicroneme, respectively 



(Fig. 2, B and C). Interestingly, when overexpressed under 
the strong constitutive GRAl promoter, a significant fraction 
of GRA24 is seen in the PV space similarly to GRAl, thus 
suggesting that GRA24 export across the PVM is a limiting 
step (Fig. 2, F and G). We concluded from these data that 
GRA24 is a novel DG-like protein that crosses the PVM of 
maturing vacuoles to reach the host cell nucleus, thereby be- 
longing to the recently defined GRA exported protein sub- 
family typified by GRA16 (Bougdour et al., 2013). 

GRA24 forms a stable dimeric complex with host p38a kinase 

Because GmV24 operates as a secreted DG that crosses the 
PVM and traffics to the host cell nucleus, we next sought 
partners of GRA24 by applying chromatography and ion- 
trap mass spectrometry to nuclear extracts that were prepared 
from inducible GRA24-HAFlag— expressing 293 cell lines. 
After induction, the protein localizes in the nucleus (Fig. 3 A), 
and regardless of the tag position (N- or C-terniinal fusion), 
GBJ\.24 is embedded in a remarkably robust molecular mass 
complex (400 kD, fractions 22—28; Fig. 3, B and E), which 
contains the CSBP2/SAPK2a isoform, also named p38a, as 



2074 



Toxoplasma GRA24 activates iiost p38 IVIAPK | Braun et al. 



Article 



A 







CM 












< 




O 
00 
3 


o 

00 











B 



Mr(k) 
39 - 

39 - 



51 



CNCNCNCN CN 



HFF 



-«*P-p38a 

-*TgHDAC3 

*P-p38a 



HFF u.i 



RHkuSO 




RHkuSO 
Agra24, gra24+ 



I, 



CM 



Mr(k) 
39 - 

39 - 

39 - 
22 - 

D 



X 

u.i. ca 
C~N c"?? 



C N C N C N 



•*P-p38a 

• p38a 

ITBP 

• Toxofilin 



J774 



< < 



3 3 3 3 3 
.S 5 - - 



X X 

en a: 



E a 8 

CL CL Z 



39 - 






64 - 




64 - 


* • 




22 - 







i*P-p38a 

*P-Ser133 
' CREB 

*P-Thr71 
ATF2 

IToxofilin 



Figure 4. GRA24 elicits phosphorylation and relo- 
cation of host p38a in the nucleus. (A and B) Immuno- 
blotting detection of *P-p38a and p38oi in HFF (A) and 
J774 M0 (B) uninfected (u.i.) or infected (24 In) witin 
KHkuSO, RHku80Agra24. PrukuSO, and Pruku80Agra24 
strains. Cytosolic (C) and nuclear (N) cell lysates were 
probed with the indicated antibodies. Histone H4 acety- 
lated (K5-K8-K12-K16) (nuclear), TBP (host-specific), and 
TgHDAC3 (parasite-specific) levels are shown as loading 
controls. (C) Subcellular in situ detection of *P-p38c( 
(green) in HFF uninfected (u.i.) or infected (24 h) with 
RHkuSO, R[iku80Agra24, RHku80Agra24, GRA24\ and 
Neospora caninum strains. Data are representative of at 
least three experiments. Bars, 10 ^m. (D) Immunoblotting 
detection of *P-p38c(, *P-Ser133-CREB, and *P-Thr71- 
ATF2 in nuclear fraction and Toxofilin (parasite-specific) 
in cytosolic fraction from HFF uninfected (u.i.) or infected 
(24 h) with the aforementioned strains and Neospora 
caninum. Data are representative of two experiments. 



the only GRA24 relevant partner (Fig. 3 C). Western blot- 
ting confirmed the partnership between GRA24 and p38a 
(Fig. 3, B and E).The GRA24— p38a complex was also iso- 
lated from primary human foreskin fibroblasts (HFFs) and 
mouse J774 M0 line hosting GRA24-HAFlag— expressing 
tachyzoites, demonstrating that the complex forms not only 
in stromal cells but also in macrophages hosting cell-cycling 
tachyzoites (Fig. 3 D). Using the dual phosphorylation of 
Thrl80/Tyrl82 as a haUmark of p38 activation (Coulthard 
et al., 2009), we confirmed by immunoblotting that p38ot in 
complex with GRA24 is activated (i.e.,phosphoThrl80/Tyrl82 
positive) when purified from inducible GRA24-HAFlag— 
expressing 293 cells (Fig. 3, B and E) or from infected murine 
and human cells (Fig. 3 D) . 

Deletion of GRA24 is sufficient to abolish 
host cell p38a activation by Toxoplasma 

To investigate the role of GRA24 on host cell p3Ba activation 
during infection, WT, GRA24 knockout, or complemented 
tachyzoites were generated.The GRA24 gene was disrupted by 



homologous recombination (Agra24) in both type I (RHfeiiSO) 
and type II (Prufew^O) strains while the mutant was comple- 
mented with a single copy of a GRA24 functional allele 
driven by its endogenous promoter (Fig. SI and Table SI). All 
strains behave similarly in vitro, indicating that QRA24 is not 
essential for either invasion or intracellular survival (unpub- 
hshed data). Consistent with our previous data, no phospho- 
p38a population is ever detected in HFF and J774 M0 cells 
hosting GRA24 knockout tachyzoites in contrast to cells 
hosting GRA24 WT virulent type I or cystogenic type II 
parasite strains (Fig. 4, A and B). Importantly, reintroduction 
of the GRA24 gene in the mutant restores a parasite's ability 
to induce p38a phosphorylation (Fig. 4, C and D), whereas 
p38a was never found to be activated by Neospora caninum re- 
gardless of the multiplicity of infection (MOI) or the duration 
of infection (Fig. 4, C and D). Remarkably, although p38a 
activation usually occurs transiently and within minutes in 
response to most stimuli (Ferreiro et al., 2010), GRA24- 
dependent p38a phosphorylation is, in contrast, impressively 
sustained as late as 24 h after the swift tachyzoite invasion step 
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(Fig. 4), a trait which fits well with the unusual long-lasting 
presence of GRA24 into the host cell (Fig. 2 A) and suggests 
a direct control on the p38a kinase activity. 

Once present in the hosting cell, GRA24 allows 
sustained p38a autophosphorylation, independently 
of the prototypic kinase cascade 

The prototypic p38a activation is achieved by upstream 
MAPK kinase (MKK3/6; Ashwell, 2006) or, alternatively, 
through TAB 1 -mediated p38a autophosphorylation (Ge 
et al., 2002). Another pathway, restricted so far to T cells 
upon T cell receptor activation, includes phosphorylation of 
Tyr323 by the ZAP-70 tyrosine kinase that leads to p38 au- 
tophosphorylation but solely on ThrlSO (Salvador et al., 
2005). The host p38a autophosphorylation, observed when 
tachyzoites are nested in their PV, was proposed to reflect a 
TAB 1 -dependent mechanism, thus a mechanism bypassing 
the MKK-dependent MAPK phosphorylation (Kim et al., 
2005). In T. ^^o/7(f/'!— hosting J774 M0 cells exposed to 
SB203580, a competitive inhibitor of ATP binding blocking 
p38a autophosphorylation but not trans-phosphorylation 
by MKKs (Gum et al., 1998), there is a complete abrogation 
of GRA24-induced p38a phosphorylation, a phenotype which 
was not observed when cells were treated with PD98059 
and SP600125, which have been shown to block the MEKl/2 
and the JNK pathways, respectively (Fig. 5, A and B). In addi- 
tion, expression of GRA24 is sufficient to trigger p38a auto- 
phosphorylation because SB203580 inhibits p38a activation 
in 293 cells induced for GRA24 synthesis, whereas MEKl/2 
and JNK inhibitors do not (Fig. 5, C and D). GRA24 is the 
first molecular switch originating from a eukaryotic parasite 
that directly activates p38a in a sustained way independendy 
of the prototypic kinase cascade in both human stromal cells 
and mouse macrophages. 

GRA24 interacts with p38a through a singular 
kinase interacting motif (KIM) 

To detail further how GRA24 binding to p38ot promotes ac- 
tivation of the kinase, we mapped the GRA24 region inter- 
acting with the enzyme. To this end, HAFlag epitope-tagged 
carboxy-terminal truncated forms of GRA24 were stably ex- 
pressed in RHfe/Y.SO parasites (Fig. 6 A). The two chimeric 
proteins, namely GRA24 AR2 and GRA24 AR1/R2, mi- 
grate at their expected molecular weights (Fig. 6 C) and lo- 
calize in the host cell nuclei after infection (Fig. 6, B and C), 
thereby indicating that the removed domains do not affect 
the export and the transit into the host cell of the chimeric 
proteins. However, although deleting only repeat R2 has no 
effect on p38a binding and activation, double deletion of Rl 
and R2 aboHshes GRA24 binding property to p38a (Fig. 6 C), 
preventing any p38a phosphorylation (Fig. 6 B), thus point- 
ing to the required contribution of GRA24 Rl domain to 
p38oi activation. These observations prompted us to look for 
the molecular determinants within the repeat that drive the 
interaction with/activation of p38a.We identified a perfect 
replicate of the canonical KIM (also known as D motif. Fig. 6 D) 
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Figure 5. GRA24 promotes p38a autophosphorylation. (A) Diagram 
showing the enzymes and their respective inhibitor. (B) Specific pharma- 
cological inhibition of GRA24-dependent p38a phosphorylation in J774 
M0. Cells were incubated for 1 h with p38 (SB203580, 15 pM), SAPK/JNK 
(SP600125, 25 pM), and MKK (PD98059, 15 pM) inhibitors or DMSG ve- 
hicle before infection in presence of the drug with RHtoSOand RHkuSO 
Agra24 strains (18 h). Nuclear fractions were immunoblotted with the 
indicated antibodies. Data are representative of two experiments. 
(C) In situ *P-p38a (green) detection in 293-TRex cells before (— Tet) or after 
(+Tet) conditional GRA24-HF expression (18 h; red). SB203580 impedes 
p38oi phosphorylation. Images are representative of at least three experi- 
ments. Bar, 2.5 pm. (D) p38a phosphorylation status in 293-TRex cells 
pretreated by inhibitors (see B) and induced (+Tet) or not (—Tet) to ex- 
press GRA24-HF. Nuclear fractions were immunoblotted with the indi- 
cated antibodies. Data are representative of two experiments. 



defined as a key MAPK-interacting structural element that, in 
addition, provides specificity toward the MAP kinase partners 
(Remenyi et al., 2005; Francis et al., 2011). The motif is typi- 
fied by a cluster of basic residues at the N terminus and hy- 
drophobic amino acids near the C terminus (Fig. 6 D) . The 
GRA24-KIM clearly belongs to the tyrosine phosphatase 
PTP-KIM class (defined in Garai et al, 2012; Fig. 6 D) but 
displays an unusual insertion of a proline (Pro-451) between 
hydrophobic canonical residues 4)L and 4)A that, in tandem 
with Pro-452, confers local rigidity to the linear peptide em- 
bedded in the KIM binding site (Fig. 6 D) and, most likely, 
adds specificity toward p38a (Garai et al., 2012). Our docking 
models vahdate the proper binding of GRA24-KIMs to p38a 
(Fig. 7). More specifically, the two N-terminal arginine resi- 
dues (Arg-444 and Arg-445) share salt bridges with CD-site 
Asp-316 and Glu-81 of p38a, whereas both Leu-450 (4>L) 
and Leu-453 (4)A) insert optimally into the shallow hydro- 
phobic pockets (Fig. 7, A and B). Finally, Ile-455 (<t)B) rein- 
forces the hydrophobic interactions between the two partners, 
as previously reported for the He-PTP Val-31 residue (4)B; 
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Figure G. A singular KIM motif in GRA24 repeats is required for 
p38a autophosphorylation. (A) Diagram of GRA24 C-terminal trunca- 
tions. Signal peptide (yellow), nuclear localization signal (green), internal 
repeats (red), and KIM (red square) sequences are shown. (B) IFA of GRA24 
chimeric proteins (red), subcellular localization, and p38a phosphoryla- 
tion status (green) in HFF infected with parasites expressing HF-tagged 
GRA24 full-length or truncates (see A). Data are representative of three 
experiments. Bars, 20 pm. (C) J774 M0 were infected (24 h) with parasites 
expressing HF-tagged GRA24 full-length or truncates. Tagged proteins 
were Flag-immunoprecipitated and eluates were analyzed by immuno- 
blotting. Data are representative of two experiments. (D) Sequence align- 
ment of KIMs from GRA24 and phosphatases. The N-terminal 
hydrophobic residue (cpH), positively charged residues (basic), and 
(pL-X-cpA-X-cpB motif are indicated. 
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Figure 7. Structural features of the GRA24-KIM1 peptide main 
chain docl<ed into p38a. (A and B) GRA24-KliVI1 peptide docked into 
p38c« (1 LEW). The N-terminal side binds to the charged CD site while 
C-terminal hydrophobic residues insert into hydrophobic pockets (pL, ipA, 
and (pB. The hydrophobic pocket is formed by the loop between P7 and 
P8 and the loop joining aD and aE helices. (Cand D) GRA24-KIIVI1 pep- 
tide docked and energy minimized (see Materials and methods) into p38 
KIM-binding domain. (C) Hydrogen bond network of the two arginines at 
the CD domain. (D) Peptide hydrophobic C terminus flanked by loops 
159-163 and 118-125. (E) Coi-backbone KIM binding site area superposi- 
tion of p38a, ERK2, and JNK models demonstrating local displacement 
upon binding. p38a peptides (MEF2A-1lew and MKK3b-1lez green), ERK2 
peptides (MKP3-2fys and HePTP-2gph, blue), and JNK (pepJIPl-lukh, 
purple) and p38a-GRA24/KIMl docked peptide (light blue) are shown. 



R2 (Fig. 6 A), might favor p38a trans-autopliosphorylation 
via the simultaneous binding of two p38a monomers (Diskin 
etal.,2007). 



Francis et al., 2011). Collectively, peculiar GRA24-KIM 
structural features empower GRA24 ability to bind to p38a, 
thereby accounting for the remarkable resistance to the highly 
stringent conditions used during the GRA24— p38a complex 
purification (Fig. 3, B and E). Such binding properties are 
expected to potentiate p38a autophosphorylation, as KIM- 
niediated interaction with p38a has been shown to trigger 
local (Fig. 7 E) as well as long-range conformational changes, 
in particular within the activation loop that becomes disor- 
dered, therefore exposing the phosphorylatable Thr-180 and 
Tyr-182 residues (Chang et al., 2002; Zhou et al., 2006; Zhang 
et al, 201 1). Of note, the KIM repeat itself, i.e., within Rl and 



GRA24 activates the expression of proinflammatory 
genes in macrophages 

The unique features of the unraveled interaction between 
GRA24 and p38a and the consequences on p38a activation 
status prompted us to analyze whether the GRA24— p38a 
complex contributes to the typical changes in gene expres- 
sion observed in either stromal or phagocytic cells infected 
with tachyzoites (Blader et al., 2001; Saeij et al., 2007; Jensen 
et al., 2011). To address this possibility, we performed a genome- 
wide expression profiling of mouse BM-derived macrophages 
(BMDMs) hosting parental versus Agra24 parasites from both 
type I and type II strains (Fig. 8 A). We focused our analysis 
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on genes that were altered with more than a twofold change 
when comparing each Agra24 mutant with their respective pa- 
rental strains. KEGG analysis highlights 10 pathways as signifi- 
candy and selectively induced (P < 0.05) in a GRA24-dependent 
fashion, a majority of which regulates proinflammatory cyto- 
kines (e.g., IL-6, IL-i2h, and TNF) and chemokines (e.g., 
MCP-1, CCL5, CXCU, and CXCLIO) known to operate m 
any sterile tissues colonized by invasive microbes, whether 
they are nonparasitic or parasitic (Fig. 8 C). Interestingly, the 
GRA24-dependent changes are exacerbated in cells infected 
with the type II clonal lineage, possibly as a consequence of 
allelic polymorphism between the strain types (Fig. 1 B), which 
leads to a concomitant up-regulation of 13 genes linked to 
M1-M0 activation and down-regulation of 17 genes asso- 
ciated with the alternative M2-M0 pathway (Fig. 8 B). 
We have also noticed that both IFN-7- and iNOS-coding 
transcript populations were also up-regulated in a GRA24- 
dependent manner by 12- and 3-fold, respectively (Fig. 8, A 
and B). Of note, at least within the time frame (24 h) ex- 
amined in these experiments, the BMDM supernatants dis- 
play no difference in IFN-7 levels as assessed by ELISA 
(unpublished data). 

GRA24-dependent IL-12p40 and MCP-1 production 
by mouse macrophages hosting tachyzoites 

As stated above, in laboratory mice, the T. gondii— and mouse- 
driven generation of immune type 1 effectors that prevent both 
uncontrolled tachyzoite population expansion and subsequent 



host death is characterized by the IL-12p40 production by 
tachyzoites-hosting macrophages (Scharton-Kersten et al., 
1996; Sher et al., 2003; Yap et al, 2006). GRA24 is clearly 
playing a major role in regulating IL-12p40 synthesis, as its 
loss severely compromises both IL-12 niRNA and protein 
levels in mouse BMDM (Fig. 8 B and Fig. 9 A). Similarly, on 
day 4 after infection, mice infected with a PruknSO Agm24 
strain had significantly less IL-12p40 in their i.p. cavities and 
serum than mice infected with the parental strain (Fig. 9 B). 
However, despite an efficient GRA24-mediated p38a activa- 
tion (Fig. 4), type I parasites fail to trigger IL-12 secretion in 
BMDM (Fig. 8 A and not depicted), a phenotype which 
might reflect a dominant IL-12 suppressive activity driven by 
the ROP16 type I allele (Saeij et al, 2007;Yaniamoto et al., 
2009; Butcher et al., 2011). In addition, other effectors, such 
as GRA15, are involved in the control of IL-12 synthesis 
through the activation of the NF-kB pathway (Rosowski et al., 
201 1). Although GRA15 does not trigger p38a phosphoryla- 
tion (Fig. 9 C), GRA24 does not affect NF-kB (not depicted), 
highlighting that several T. gondii strain-specific or -shared 
effectors display synergistic or antagonist effects on the regu- 
lation of IL-12p40— coding gene expression. In support of 
GRA24 key contribution to the IL-12 response, GRA24 is 
found to critically account for the typical increase in MCP-1, 
CXCLl, CCL5, and CXCLIO chemokine synthesis (Fig. 9 D) 
triggered by a type II tachyzoite lineage. Remarkably, MCP-1 
was shown to play a substantial role in the recruitment, to 
T. gondii— loaded tissues, of Gr-1^ inflamniatory monocytes 
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Figure 8. GRA24 alters the host-cell transcrlp- 
tome. C57BL/6 BMDMs were infected with either type I 
RH/ru80Agra24or type II PrukuSO Agra24 mutant para- 
sites and compared with infection with their respective 
parental Toxop/osmo strains. Cells were infected at an 
MCI of 1:3 and total RNAswere extracted 18-20 h after 
infection. cDNA synthesis, labeling, and differential gene 
expression analysis are detailed in Materials and meth- 
ods. The differentially regulated genes were identified 
after normalization of the raw data and filtering genes 
with a fold change cutoff of twofold or more and p-value 
<0.05. (A) Heat map of the GRA24-regulated cytokine- 
related genes identified by KEGG analysis. Mean logj 
gene expression values were median-centered and genes 
were clustered according to the biological pathways 
identified. The complete set of genes is listed in GEO 
dataset accession no. GSE38782. (B) Heat map represen- 
tation of the GRA24-regulated genes associated with 
classical (Ml) or alternative (M2) M0 phenotypes after 
infection of BMDM with type II strains. Gene symbols 
and the fold change (EC, at least twofold mean differ- 
ence) in gene expression are listed for GRA24ll-depen- 
dent genes for comparison with those reported by Jensen 
etal. (2011) as R0P16l-and GRAl5ll-regulated. (C) KEGG 
analysis of the differentially expressed genes when com- 
paring WT versus Agira24-infected cells. Genes listed in 
GEO dataset accession no. GSE38782 (fold change cutoff 
twofold or more) were analyzed to determine pathways 
that were statistically overrepresented (p-value <0.05). 
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Figure 9. GRA24 promotes chemokine secretion and con- 
tributes to tPie control of early parasite replication in vivo 
at the site of infection. (A) IL-12p40 cytokine production 
by uninfected (u.i.) or 24-h-infected {Pruku80Agra24 or paren- 
tal strains) BMDM measured using ELISA. Means of tinree inde- 
pendent experiments ± SD are sinown (*, P < 0.005, two sample 
Student's f tests). (B) Peritoneal lavage fluid and serum were 
collected on days 2 and 4 after infection of C57BL/6 mice that 
had received an i. p. dose of 10^ Pru/tuSOWT or PrukuSO Agra24 
tachyzoites. Concentrations of IL-12p40 were determined by 
ELISA. Data shown are means ± SD with n = 3 individual mice 
per parasite genotype at each time point. Error bars represent SD 
from one experiment. *: significantly different from WT (P < 
0.05, Student's ttest). (C) Subcellular detection in J774 infected 
with Agra24ar\d 4gra 75 mutants demonstrated that *P-p38a 
and Egr-1 induction is dependent of GRA24 and independent of 
GRA15. Data are representative of two experiments. (D) Com- 
parison of GRA24-dependent and -independent chemokine ex- 
pression profiles after BMDM infection by Toxoplasma. 
Quantification is given in digital light units (DLU) drawn from 
the whole chemokine array. Error bars represent SD (two sample 
Student's t tests; *, P < 0.006). (E) WT C57BL/6 mice were i.p. 
infected with a dose of 10^ Pru<ru80WTor PrukuSO \gra24 
tachyzoites. Parasites were enumerated using recovered i.p. con- 
tents or from homogenates of spleen at the indicated times 
after infection. Data shown are means + SD (from three mice 
per parasite genotype) and representative of two independent 
experiments. *: significantly different from WT (P < 0.05, Stu- 
dent's ttest). 



tliat, in turn, produce nitric oxide and are essential for the 
local control of the tachyzoite population size (Robben et al., 
2005; Dunay et al., 2008). The above results argue that the 
GRA24-dependent MCP-1 secretion might be involved in 
the recruitment of Gr-1* monocytes in tissues transiently 
remodeled as cell-cycling tachyzoite-protective niches. To 
examine the functional consequences in vivo of GRA24- 
dependent MCP-1 secretion, we challenged C57BL/6 mice 
with 10^ tachyzoites of FrukuSOWT or PrufeK^O Agra24 
strains and assessed parasite numbers in the peritoneal exu- 
date as well as in the spleen. A larger number of parasites are 
detected beyond day 7 in the peritoneal contents of mice 
infected with PrukuSO Agra24 than in mice infected with 
the parental strain (Fig. 9 E). Parasites are also found in the 
spleen, indicating that the infection becomes quickly sys- 
temic, but in this tissue no significant deviation in parasite 
loads was detected between the WT and the Agra24 strains at 
any of the time points analyzed (Fig. 9 E). These results sug- 
gest that the deletion of GRA24, while suppressing MCP-1 
production, most likely prevents Gr-1"'" recruitment and con- 
tributes to impair the ability of the mice to suppress parasite 
numbers at the site of infection in the gut. We then assess if 
GRA24 is altering the overall host survival of C57BL/ 6 mice. 
A dose of 5 x 10^ tachyzoites required to achieve a near LDjqo 
when using the Fiuku 80 strain was i.p. inoculated to C57BL/6 



mice whose survival was monitored over the acute phase of 
infection. Mice infected by either the WT Prufeti^O or VrukuSO 
Agra24 parasites succumb to infection within 10 to 20 d (un- 
published data; n = 6 mice per strain). Therefore, QRA24- 
deficient parasites did not fully recapitulate the phenotype of 
MCP-1~^^ and CCR2^^^ mice that were shown to be pro- 
foundly susceptible to acute toxoplasmosis (Robben et al., 2005; 
Dunay et al., 2008). This difference might be a result of the 
genetic background of the parasite strains, the type II ViukuSO 
strain (Fox et al., 2011) being more virulent and less cysto- 
genic than the type II ME49 strain used in most Toxoplasma 
infection including those using MCP-1 and CCR2^^^ mice 
(Robben et al.,2005). 

GRA24-p38a complex regulates the expression 
of Egr-1 and c-Fos transcription factors 

Finally, as a link between the cytokine response and the direct 
GRA24-dependent p38a activation, we identified several 
transcription factors up-regulated in a GRA24-dependent 
fashion in M0, among which those that either control the 
immediate early gene response or are immediate early proto- 
type gene products hke Egr-1 and c-Fos (Fig. 10 A). Of note, 
first, Egr- 1 and c-Fos belong to a restricted set of 30 core 
genes whose expression is equally regulated by major p38a 
MAPK-activating stimuli (i.e., osmo-stress, the cytokine TNF, 
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Figure 10. GRA24 triggers up-regulated 
expression of the host transcription fac- 
tors Egr-1 and c-FOS. (A) Heatmap repre- 
sentation of host transcription factors 
up-regulated in a GRA24-dependent fasliion 
in BMDIVl infected by tine aforementioned 
strains. Egr-1 and c-FOS are indicated in red. 
(B) IFAof GRA24-dependent induction of Egr-1 
(red) in uninfected (u.i.) or infected (24 h) 
HFF by the indicated strains. Images are rep- 
resentative of at least four experiments. Bar, 
20 pm. (C) Cellular fractionation of J774 cells 
left uninfected (u.i.) or infected (24 h) by the 
indicated strains. Cytoplasmic (C) and nuclear 
(N) fractions were analyzed by Western blot- 
ting using indicated antibodies. TBP (host- 
specific) and toxofilin (parasite-specific) are 
shown as loading controls. Data are represen- 
tative of three experiments. (D) Western blot 
analysis of nuclear fractions of 293-TRex cell 
line expressing 6RA24-HF using indicated 
antibodies. Cells were either left untreated 
(— Tet) or treated (+ Tet) with 1 pg/ml tetra- 
cycline for 20 h. AcH4 was used to assess host 
nuclear protein enrichment. Data are repre- 
sentative of three experiments. (E) IFA show- 
ing the serine 32 phosphorylation status of 
c-FOS (in red) in serum-starved confluent HFF 
left uninfected (u.i.) or after 24 h of infection 
by RHkuSO, RHkuSO Agra24. and RHkuSO 
Agra24, GRA24+ parasites. Images are repre- 
sentative of two experiments. Bar, 20 pm. 



and the protein synthesis inhibitor anisomycin; Ferreiro et al., 
2010), and second, the Egr-1 promoter activity is blocked by a 
pharmacologically inactivated p38a kinase (Tur et al., 2010). 
These observations strongly argue for a direct function of 
GRA24 on classical p38a MAPK-regulated genes and con- 
firm previous studies on the induction of Egr-1 and its paralog 
Egr-2 during Toxoplasma infection (Phelps et al., 2008; Wiley 
et al., 2011). By Western blotting and immunofluorescence, 
we indeed confirmed at the protein level that Egr-1 is induced 
in the nuclei of human stromal cells and mouse macrophages 
hosting WT T. gondii, a phenotype which is not induced in 
cells hosting Agra24 Toxoplasma or Neospora (Fig. 10, B and C; 
and not depicted). Reintroduction of the GRA24 gene re- 
stored the Toxoplasma ability to up-regulate Egr-1 (Fig. 10 B). 
Egr-1 induction was also vahdated in 293 cells conditionally 
expressing GRA24 (Fig. 10 D), and the level of expression 
was reduced in a dose-dependent fashion by the inhibitor 
SB203580, thereby establishing the link between Egr-1 and 
the autophosphorylation of p38a by GRA24 (Fig. 5 B). Fi- 
nally, we identified the GRA24-KIM as the basis of p38a 
binding and autophosphorylation that further entails Egr-1 
activation (Fig. 6 B). Besides Egr-1, c-fos is also regulated at the 
transcriptional level by GRA24 (Fig. 10 A). In particular, this 
transcription factor is phosphorylated on serine 32 in a GRA24- 
dependent fashion in both murine M0 and human fibroblasts 
(Fig. 10 E and not depicted). Unexpectedly, Toxoplasma elicits 



a prolonged induction of Egr-1 and c-Fos, which correlates 
with the nuclear persistence of the GRA24— p38ot complex 
during infection but differs from the transient usual canonical 
activation of these genes. The induction of these transcription 
factors by GRA24 is most likely followed by a second wave 
of host cell transcriptome alteration because both Egr-1 and 
c-Fos then operate in a hierarchical manner by inducing the 
expression of their downstream genes (e.g., Egr-1 -regulated 
genes, see GEO accession no. GSE38782). 

DISCUSSION 

Host cell invasion is a critical step by which Toxoplasma zoites 
rapidly establish themselves and further complete their devel- 
opmental program as bradyzoite-containing cysts known to 
be located in the tissue of prey that are ingested by predators. 
Recent research has highlighted that early during the inva- 
sion process, the parasite secretes several rhoptry bulb resident 
proteins with a majority of kinases and phosphatases that co- 
opt host cells by interfacing with their signaling pathways 
(Melo et al., 2011). Once internalized, a Toxoplasma tachyzo- 
ite resides within a dynamic PV that is constantly reshaped 
through molecular exchanges with the host cell cytoplasm 
and that support parasite growth. Although nonenzymatic 
molecules secreted by DG were first thought of as parasite 
effectors restricted to the PV biogenesis and maturation, a re- 
cent study (Bougdour et al., 2013) and the present study 
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highlight that DG-resident nonenzymatic proteins are se- 
creted into the vacuole space, cross the PVM, and reach the 
host cell nucleus where they modulate host gene expression. 
Thus, GRA16 and GRA24 emerge as a subfamily of DG- 
scaffolding proteins exported beyond the tachyzoite-hosting 
vacuole to operate as effector molecules in infected cells. 

In the present study, we have identified GRA24 as a strong 
interactor of host p38a, with the resulting complex acting as 
a regulator of gene expression in the host cell nucleus. Indeed, 
we establish that GRA24 interacts with the p38a enzyme, 
triggers its autophosphorylation, and consequently its activa- 
tion. GRA24 bypasses the typical MKK-dependent MAPK 
phosphorylation cascade upstream of p38 and is the first parasitic 
microbe agonist so far described to promote p38a autophos- 
phorylation through a dimeric stable complex. Interestingly, 
GRA24 specifically promotes ThrlSO autophosphorylation 
of p38a because it neither activates ERK1,2 (Thr202/Tyr204 
and Thrl85/Tyrl87, respectively) nor alters theJNK (Thrl83/ 
Tyrl85) phosphorylation (unpublished data). In agreement 
with these results, none of the aforementioned kinases copu- 
rify with GRA24, reinforcing the notion of a high specificity 
of GRA24 function on p38a. Another unique feature of 
GRA24 function lies in its property to promote a sustained 
phosphorylation of p38a that correlates with the long-lasting 
secretion of GRA24.This kinetics of activation is remarkably 
distinct from the transient activation mediated by the proto- 
typic MAPK pathway and has important consequences be- 
cause it likely avoids conventional negative feedback by 
phosphatases. Consistent with this idea, we provide evidence 
that the carboxy-terminal repeat of GRA24 is essential for 
binding to p38a and its subsequent autophosphorylation 
(Fig. 6 B). Furthermore, molecular modeling reveals that 
GRA24 repeat contains a p38a docking KIM motif that is 
phylogeneticaUy related to those of the PTP tyrosine phos- 
phatase family (Fig. 6 D). Collectively, these data suggest that 
GRA24 has evolved a structurally optimized KIM domain 
that might be able to compete accurately with the phospha- 
tases for binding to MAPK. 

Toxoplasma-induced p38a activation was initially reported 
as TAKl- and Myd88-iiidependent and proposed to occur 
through an association with the scaffold protein TAB 1 (Kim 
et al., 2005). But — of note — first, TAB 1 does not activate p38a 
after TNF exposure of fibroblasts derived from the MKK3- 
MKK6 double knockout mice (Brancho et al., 2003), and 
second, once bound to p38a,TABl prevents the p38a nu- 
clear localization (Lu et al., 2006) disproving the validity of 
this proposal. In contrast, our many complementary data allow 
us to highlight a major role for GRA24 as the shuttling and 
scaffolding protein that allows not only the p38a nuclear lo- 
calization but also its autophosphorylation, at least in two dis- 
tinct cell Uneages — fibroblasts and macrophages — that T. gondii 
tachyzoites do subvert as host cells. Moreover, a contribution 
of parasite-derived tyrosine kinase is also very unlikely be- 
cause bypassing the recourse to tachyzoites, through the ectopic 
expression of GRA24 in mammalian cells, is clearly sufficient 
to promote p38a autophosphorylation. It is noteworthy that 



GRA24 has multiple alternative splice variants. Specifically, the 
presence of the exon E3b in the niRNA isoform 2 introduces 
a premature stop codoii, leading to a protein variant devoid of 
the terminal repeats Rl and R2 (Fig. 1 A). In this fashion, the 
parasite might use alternative splicing to turn off GRA24- 
dependent p38a autophosphorylation under certain growth or 
developmental conditions, which remain to be discovered. 

The ability of GRA24 to directly promote host p38a 
autophosphorylation places it as a potential regulator of the 
synthesis of growth factors and/ or cytokines/ chemokines 
and associated receptors, at least in macrophages once they are 
subverted as tachyzoite-hosting cells. Using microarray analy- 
sis, we showed that GRA24 controls, at the transcriptional 
level, an impressive network of genes encoding for cytokines/ 
chemokines (Fig. 8), well beyond those associated previously 
with the secretion of GRA15 (Jensen et al., 2011; Rosowski 
et al., 2011). Chief among the cytokines is the proinflamma- 
tory cytokine IL-12p40, with a 20-fold transcript excess that 
typifies a robust Ml profile in M0 hosting type II tachyzo- 
ites. As reported previously, IL-12 synthesis induction by 
Toxoplasma-inkcted MO occurs at lower levels and with dis- 
tinct kinetics relative to endotoxin stimulation; in brief, al- 
though IL-12 levels peaks within 6 h of LPS stimulation (Kim 
et al., 2005), parasite-induced IL-12 is detected within a 24— 
40-h window after macrophage invasion (Fig. 9 A and Kim 
et al., 2005). This timing is consistent with the kinetics of 
GRA24 secretion and the requisite for de novo synthesis of 
p38-activated transcription factors driving IL-12 synthesis. Of 
note, once macrophages are hosting tachyzoites, both IL-1(3 
and IL-18 signaling use MyD88 as an adaptor, leading to 
NF-kB and MAPK activation and subsequent IL-12 produc- 
tion (Adachi et al., 1998). It is also noteworthy that IL-18 is 
approximately sixfold induced by GRA24 in MO infected 
with type II strain (Fig. 8 A). It is therefore possible that auto- 
crine production of IL-18 results in GRA24-mediated IL-12 
release. Anyhow, the massive production of IL-12 induced by 
type II strains when tachyzoites start to proliferate is most 
likely controlled by the synergistic action of GRA15 and 
GRA24, which keep parasite numbers low, thus facilitating 
long-term parasitism. This is consistent with the observation 
that GRA24 is one among the few polymorphic highly ex- 
pressed genes that are under positive selection in Toxoplasma 
population (Minot et al., 2012). We can also speculate that 
GRA24 loss in Neospora caninum correlates with the inability 
of Neospora to properly modulate the host immune response, 
an argument which has been put forward earlier to explain 
why Neospora has a more hniited host range than Toxoplasma 
(Reid etal.,2012). 

Besides triggering a strong Thl response, GRA24 also sig- 
nificantly contributes to suppress M2 markers at the tran- 
scriptional level in MO hosting type II tachyzoites, a situation 
which powers the Ml response. The repertoire of M2 genes 
that are down-regulated in a GRA24-dependent manner 
again significantly exceeds that of GRA15, with Argl being 
distinctively ^60-fold repressed by GRA24 Qensen et al., 
2011; Fig. 8 B). However, neither arginase-1 nor the mannose 
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receptor type C (Mrcl/CD206) are induced at the protein 
level in A_^ra24-infected BMDMs (unpublished data), sug- 
gesting that the suppression of the GRA24-dependent tran- 
scriptional repression of M2 genes in type II strain is not 
sufficient to induce the production of the corresponding pro- 
tein product, consistent with previous studies (Butcher et al., 
2011;Jensen et al, 2011). Additionally, GRA24 is found to 
critically account for the typical increase at the protein level 
of four cheniokines including MCP-1/CCL2 (Fig. 8 A and 
Fig. 9 D). Remarkably, in niice lacking MCP-1 or its receptor 
CCR2, no recruitment, in T. gondii— hosting tissues, of inflam- 
matory Gr-l"*" monocytes, is detectable, a situation which Ukely 
accounts for the sustained parasite burden and early host 
mortality (Robben et al., 2005; Dunay et al., 2008). Consis- 
tent with these data, we show that mice infected with Agra24 
tachyzoites fail to control parasite numbers in the gut (Fig. 9 E), 
reveahng a contributory role for GRA24 in resistance to acute 
toxoplasmosis, most likely through the regulation of cytokine 
expression and secretion including MCP-1. 

GRA24 also significantly induces the secretion of CXCLIO/ 
IP- 10 (Fig. 8 A and Fig. 9 D).This chemokiiie is required for 
preventing uncontrolled T. gondii tachyzoite burden (Khan 
et al., 2000), and its depletion prevents T cell recruitment and 
effector function in mouse tissues where the tachyzoite to 
bradyzoite developmental transition is known to occur. In 
this regard, it was recently shown that the CXCLlO/IP-10 
chemokine enhances the ability of CD8^ T cells to control 
tachyzoite reexpansion in cyst-hosting mouse brains (Harris 
et al., 2012). Collectively, these studies underscore the impor- 
tance of GRA24 in orchestrating a cell response during Toxo- 
plasma infections by promoting the induction of a selective set 
of chemokines, which contributes to the recruitment of in- 
flammatory monocytes previously described as a first line of 
defense controlling intestinal pathogens, including T. gondii 
(Dunay etal, 2008). 

Given that GRA24 was involved in nonprototypical di- 
rect p38a activation, we next sought to trace the factors that 
act downstream of this unusual signaling pathway. This led us 
to show that several transcription factors, including Egr-1 and 
c-Fos (Fig. 10 A), inhabit key downstream positions within 
this coordinated acute inflammatory process driven by tachy- 
zoites nested within macrophage PV Several factors are needed 
for full, regulated transactivation of the Egrl (Tur et al., 2010) 
or c-Fos (O'Donnell et al., 2012) genes. As such, the mitogen- 
and stress-induced phosphorylation of CREB (CRE binding 
protein) at Serl33, especially by p38a, has been linked to the 
transcription of several inmiediate early genes, including c-fos, 
junB, and Egr-1 (Tur et al, 2010; O'Donnell et al, 2012). In 
this regard, the Toxoplasma tachyzoite, in a GRA24-dependent 
fashion, positively affects the transcription of two CREB or- 
thologs, I.e., CREB2/ATF2 and CREBL2 (Fig. 10 A), as well 
as CREB-dependent downstream genes (GEO accession no. 
GSE38782). We first thought that GRA24 might activate 
p38a, which in turn directly phosphorylates CREB-related 
factors in macrophages hosting tachyzoites. But this is appar- 
ently more complicated because Toxoplasma induces the 
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phosphorylation of both CREB-Serl33 and ATF2-Thr-71 
independently of GRA24 (Fig. 4 D). Therefore, GRA24 
might act upstream in the pathway. It was previously shown 
that the phosphorylation of Ser 133 facilitates the ability of 
CREB to interact with CBP (CREB-binding protein) and its 
paralog p300, therefore leading to the transcription of CREB- 
responsive genes (Johannessen et al., 2004). CBP and p300 
possess an intrinsic acetylation activity and are thought to 
acetylate specific histone tails, thereby assuring relaxation of 
the chromatin environment Qohannesseii et al., 2004). As 
such, the sustained presence of GI\A24-p38a in the nucleus 
of the tachyzoites-hosting macrophages might reflect the 
need for the complex to be embedded in the chromatin to 
operate in the vicinity of p38a MAPK-responsive genes. 
A GRA24-dependent response might alternatively be medi- 
ated through Egr-1 that physically and functionally interacts 
with CBP/p300 to modulate gene transcription (Silverman 
et al., 1998). Although interesting, this hypothesis deserves 
further investigation. 

In conclusion, GRA24 extends the scope of the function 
of T. gondii DG proteins beyond their initial and dedicated 
role in PV formation and maturation. Moreover, the specific 
attributes of the GRA24-dependent p38a activation unveil a 
novel regulatory path for p38a-controlled transcriptional 
changes. This path not only shortcuts the canonical MAPK 
signaling cascade turned on in response to various stress- 
related stimuli and intracellular pathogens (Ferreiro et al., 
2010; Krachler et al., 2011) but it also differs from the alter- 
nate p38 activation pathways described so far involving the 
scaffolding protein TABl or the T cell— restricted kinase 
ZAP-70 (Brancho et al, 2003; Salvador et al, 2005). 

MATERIALS AND METHODS 

Parasites and host cells. HFF primary cells, 293~TRex, and J774 cell lines 
were cultured in DMEM (Invitrogen) supplemented with 10% heat inacti- 
vated FBS (Invitrogen), 10 niM Hepes buffer, pH 7.2, 2 niM L-glutamine, 
and 50 |.ig/ml peniciUin and streptomycin (Invitrogen). Cells were incubated 
at 37°C in 5% CO2. BMDMs were obtained from female C57BL/6 mice. 
Bone inarrow was isolated by flushing hind tibias and femurs using a 25-gauge 
needle followed by passages through an 18-gauge needle to disperse cell 
clamps. Cells were suspended in DMEM supplemented with 10% heat- 
inactivated FBS and 20 ng/ml recombinant M-CSF (Invitrogen), and incu- 
bated in a tissue culture-treated flask for 8—12 h. Then, nonadherent cells 
were harvested and transferred in 55 cm^ non— tissue culture— treated plates 
(Corning) at 4—6 X 10^ cells per plate and further incubated at 37°C, with 
5% CO2 in humidified air. After 6 d, cells were washed with PBS to remove 
nonadherent cells, harvested by dislodging with a cell scraper in ice-cold 
PBS, and replated for the assay. This method yielded a highly pure population 
of F4/80"'' macrophages by IFA.The Toxoplasma strains used in this study are 
listed in Table SI. All T. gondii strains were maintained in vitro by serial pas- 
sage on monolayers of HFFs. 

Reagents. Antibodies against HA (3F10; Roche), MIC2 (provided by 
D Sibley, Wasliington University School of Medicine, St. Louis, MO),Toxofilin 
(provided by I.Tardieux, INSERM, U1016, Institut Cochin, Paris, France), 
GRAl and GRA7 (provided by J.-E Dubremetz, UMR 5235 Centre 
National de la Recherche Scientrfique, MontpeUier, France), p38 MAPK 
(Cell Signahng Technology), Phospho-p38 MAPK (Cell Signahng Technol- 
ogy), Phospho-c-Fos (Ser32; Cell Signaling Technology), H4 Acetylated 
(EMD Millipore),TgHDAC3 (Bougdour et al.,2009),EGRl (Cell Signaling 
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Technology), and anti-TBP (Abeam) were used in the immunofluorescence 
assay and/or in Western blotting. Immunofluorescence secondary antibodies 
were coupled with Alexa Fluor 488 or Alexa Fluor 594 (Invitrogen) . Second^ 
ary antibodies used in Western blotting were conjugated to alkaline phospha- 
tase (Promega).The inhibitors SB203580, SP600125, and PD98059 were 
purchased firom InvivoGen. 

Immunofluorescence microscopy. Immunofluorescence assays were per- 
formed as described in Bougdour et al. {2009). Cells grown on coversHps 
were washed in PBS and fixed for 20 min at room temperature with PBS 
containing 3% (vol/vol) formaldehyde. Fixed cells were permeabilized with 
PBS-0.1% Triton X-100 (vol/vol) for 10 min and blocked in PBS-3% BSA 
for 1 h at room temperature. Samples were incubated in PBS containing 
3% BSA with the primary antibodies indicated in the figures, followed by 
the secondary antibodies coupled with Alexa Fluor 488 or Alexa Fluor 568 
(Invitrogen) at a 1:1,000 dilution each in PBS~3% BSA. Nuclei of both host 
cells and parasites were stained for 10 min at room temperature with Hoechst 
33258 at 2 |j,g/ml in PBS.After four washes in PBS, coverslips were mounted 
on a glass shde with Mowiol mounting medium. Images were acquired with 
a fluorescence microscope (Axioplan 2; Carl Zeiss). 

Plasmid constructs. To construct the vector pLIC-GR/124-HF, the cod- 
ing sequence of GRA24 was amplified using primers LICF-117710_F and 
LICR-117710_R using RH genomic DNA as template. The resulting 
PCR product was cloned into the pLlC-HF-dhfr vector (Bougdour et al., 
2013) using the LIC cloning method previously reported (Huynh and 
Carruthers, 2009). The C-terminal truncated versions of GRA24 were 
amplified and cloned in the aforementioned pLlC-HF-dhfr vector to gen- 
erate pLIC-GiM24AR2-HF and pLIC-GiM24ARl/R2-HEThe pLIC- 
PGRA24"Gi^24-HF was obtained by amplifying the promoter and coding 
region of GRA24 (strain GTl, Chromosome VIII, position 664,928— 
671,201) with LICF-5' region-1 17710_F2 and LICR-117710_R2 prim- 
ers. The resulting PCR product was cloned in the aforementioned pLIC- 
HF-dhfr vector. To construct the mammalian expression vector pcDNA- 
Gi?^24-HF, the coding sequence of GRA24 (residues 25—542) was 
amplified from a poly(A)"^ enriched RH cDNA library using pcDNA4- 
LIC-117710-HAFla^F and pcDNA4-LIC-l 17710-HAFla^R primers. 
The resulting PCR product was cloned by LIC into the pcDNA-LIC-HF 
plasmid (Bougdour et al. ,2013). Likewise, pcDNA4-FlagHA-LIC-l 1 77 1 0_F 
and pcDNA4-FlagHA-LIC-117710_R primers were used to generate the 
plasmid pcDNA-HF-Gi^24 harboring an N-terminal HAFlag tag. The 
vector pDEST14 KOGi^24 was generated to construct the deletion/in- 
sertion mutation of GRA24 in Toxoplasma strains RH^hSO (type I, Huynh 
and Carruthers, 2009) and PrufewSO (type II, Fox et al., 201 l).The Multisite 
Gateway Pro 3-fragment Recombination system was used to clone the dhfr 
cassette flanked by the 5' and 3' surrounding regions of GRA24 coding 
sequence of RH genomic DNA. The 5' flanking region o£ GRA24 of RH 
strain was amplified using primers attBl-117710_F and attB4-l 17710_R, 
and was cloned into the plasmid pDONR221 P1-P4 (Invitrogen). The 
3' flanking region of GRA24 was ampHfied using primers attB3-117710_F 
and attB2-117710_R, and was cloned into the plasmid pDONR221 P3-P2 
(Invitrogen). The resulting vectors, pDONR221/5'GiM2'/ and pDONR221/ 
3'Gi^424, respectively, were then recombined with pDONR221/DHFR 
(Bougdour et al., 2013) into the destination vector pDEST14, yielding 
the pDEST14 KOGiM24. The vector pDEST14 KOGiM?5 was gener- 
ated to construct the deletion/insertion mutation of gralS in T. gondii 
Prufc/(SO strain. The 5' flanking region of GRA15 of type II Pru strain was am- 
plified using primers attBl-275470_F and attB4-275470_R, and was cloned 
into the plasmid pDONR221 P1-P4 (Invitrogen). The 3' flanking region of 
gralS was amplified using primers attB3-275470_F and attB2-275470_R, 
and was cloned into the plasmid pDONR221 P3-P2 (Invitrogen). The result- 
ing vectors, pDONR221/5'^ra 15 and pDONR221/3'^rai5, respectively, were 
then recombined with the pDONR221/DHFR (Bougdour et al., 2013) 
into the destination vector pDEST14, yielding the pDEST14 KOGRAIS. 



Toxoplasma transfection. Vectors were transfected into RHfc((50 and 
VrukuSO parasites by electroporation. Electrop oration was done in a 2-nim 
cuvette in a BTX ECM 630 (Harvard Apparatus) at 1,100 V, 25 n, and 
25 |xF. Stable integrants were selected in media with 1 p.M pyrimethamine 
and cloned by limiting dilution. 

Endogenous gene tagging in Toxoplasma. We used the RHfewSO strain 
(Huynh and Carruthers, 2009) to endogenously tag GRA24 in Toxoplasma. 
In brief, ^^20 [jg of each DNA construct was Hnearized within the region of 
homology. Linearized vectors pLIC-GRy424-HF (BsiWI), pLIC-GiM24 
AR2-HF (BlpI), and pLIC-GiM24ARl/R2-HF (BlpI) were phenol-chlo- 
roform extracted and ethanol precipitated. Constructions were resuspended 
inTlowE buffer (10 mMTris-HCl, pH 8.0, and 0.1 iiiM EDTA) for transfec- 
tion. Stable recombinants were selected with pyrimethamine, single-cloned 
by limiting dilution, and verified by IFA.The resulting strains are RHfewSO 
GRA24-H¥,KHku80 GiM24AR2-HF, and KHkuSO CRyl24ARl/R2-HF, 
respectively (Table SI). 

Generation of knockout and complementation. To construct the 
deletion/insertion mutation of GRA24 in type I (RH/?/iSO) and type II 
(FrukuSO) strains of Toxoplasma, the plasmid pDEST14 KOGi^24 was am- 
phfied by PCR using primers attBl-117710_F and attB2-117710_R. After 
sodium acetate/ethanol precipitation, DNA was resuspended inTlowE buf- 
fer (10 niM Tris-HCl, pH 8.0, and 0.1 niM EDTA) and ^20 [ig of PCR 
product was used for transfection. Recombinants were single-cell cloned by 
limiting dilution and verified by PCR analysis as follows. PCR with either 
a forward primer (5' region— 1 1771 0_F2) upstream of the 5' flanking region 
used to generate the knockout construction and a reverse primer within the 
dhfr cassette (dhfr-KO-oll_R) or a forward primer within the dhfr cassette 
(dhfr-KO-oll_F) and a reverse primer (3' region— 117710_R1) downstream 
of the 3' flanking region of GRA24 confirmed a disruption in the GRA24 
locus. Further PCRs were performed to confirm the inability to amplify 
GRA24 coding sequence using primers (pcDNA4-FlagHA-LIC-117710_F 
and pcDNA4-FlagHA-LIC-117710_R) in both type I and type II strains. 
The deletion/insertion mutation of GRA15 in VrukuSO was constructed as 
described above with the plasmid pDEST14 K.OGRA15 amplified by PCR 
using primers attBl-275470_F and attB2-275470_R. Disruption in the 
GRA15 locus was confirmed by PCR analysis with either a forward primer 
(5' region— 275470_F) upstream of the 5' flanking region used to generate 
the knockout construction and a reverse primer within the dhfr cassette 
(dhfr-KO-oll_R) or a forward primer within the dhfr cassette (dhfr-KO- 
oll_F) and a reverse primer (3' region— 275470_R) downstream of the 
3' flanking region o£ gral5.To generate a complemented KMkii80Agra24 
strain, vector pLIC-PGRA24-GRy424-HF was co-electroporated with the 
pMiniHX at 1:10 ratio, followed by selection with mycophenolic acid and 
xanthine, and cloned by limiting dilution. The resulting strain is named 
KHktiSO Agra24, GRA24^ (Table SI). 

293-TRex transfection. 24 h before transfection, the cells were plated 
(80% confluency) in 6-well tissue culture dishes. 1 |j,g Flag— fusion protein- 
expressing plasmids and 0.1 )j.g puromycin selection plasmid were co- 
transfected into 293-TRex cells with Lipofectaiiiine reagent (Invitrogen) 
according to the manufacturer's instructions. 72 h later, cells were diluted in 
the presence of 5 ^ig/nil puromycin ( Sigma- Aldrich) for selection. Individual 
drug-resistant clones were expanded and tested for tetracycline-inducible 
gene expression. 

Chromatographic purification of GRA24-containing complexes. 

Nuclear extracts from 293-TRex cells stably expressing Flag-tagged protein 
were incubated with anti-FLAG M2 affinity gel (Sigma -Aldrich) for 1 h at 
4°C. Beads were washed with 10-column volumes of BC500 buffer (20 niM 
Tris-HCl, pH 8.0, 500 niM KCl, 10% glycerol, 1 niM EDTA, 1 niM DTT, 0. 1% 
NP-40, 0.5 niM PMSF, and 1 i^g ml~^ each of aprotinin, leupeptin, and pep- 
statin). Bound peptides were eluted stepwise with 250 |-Lg ml~^ FLAG peptide 
(Sigma-Aldrich) diluted in BClOO buffer. For size-exclusion chromatography, 
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protein eluates were loaded onto a Superose 6 HR 10/30 column equili- 
brated with BC500. Flow rate was fixed at 0.35 ml/niin, and 0.5-nil fractions 
were collected. 

Mass spectrometry and peptide sequencing. Protein bands were ex- 
cised from colloidal blue— stained gels (Invitrogen), treated with DTT and 
iodoacetamide to alkylate the cysteines, and then immediately subjected 
to in-gel tryptic digestion. Peptides were extracted with 5% vol/vol formic 
acid solution and acetonitrile, and injected into an Ultimate 3000 (Dionex) 
nanoLC system that was directly coupled to a LTQ-Orbitrap mass spectrom- 
eter (Thermo Fisher Scientific). MS and MS/MS data were acquired using 
Xcalibur (Thermo Fischer Scientific) and processed automatically using 
Mascot Daemon software (Matrix Science). Tandem mass spectra were 
searched against a compiled T. goudii database using the MASCOT program 
(Matrix Sciences). 

Cell fractionation. For cytosolic and nuclear fractionation analysis, cells 
were washed in ice-cold PBS and harvested using a cell scraper. After centrif- 
ugation, cells were resuspended in buffer D (10 niM Hepes, pH 7.9, 10 niM 
KCl, 1.5 mM MgClz, 0.34 M Sucrose, 10% Glycerol, 1 niM DTT, and Ix 
Roche protease inhibitor cocktail) and lysed for 8 min with Triton X-100 at 
a final concentration of 0.1%. Cell lysates were centrifuged at 1,300 g for 
5 min and the supernatant containing the cytosohc fraction was collected. 
The pellet, containing the nuclear fraction, was washed once in buffer D and 
resuspended in cell extraction buffer (Invitrogen) for protein extraction. The 
supernatant and pellet fractions were clarified by centrifugation at 20,000 g 
for 10 min at 4°C. Supernatants were mixed with protein sample buffer 
(Invitrogen) for Western blot analysis. 

Western blot. Proteins were separated by SDS-PAGE, transferred to a poly- 
vinylidene fluoride membrane (Immobilon-P; Millipore) by hquid transfer 
and Western blots were probed using appropriate primary antibodies fol- 
lowed by phosphatase-conjugated goat secondary antibodies (Promega). 
Signals were detected using NBT-BCIP (Amresco). 

In vitro cytokine ELISA. C57BL/6 BMDMs were seeded (4 x 10^ per 
well) in 96-well plates and left to adhere overnight at 37°C in 5% CO2. Cells 
were infected with freshly lysed T. gondii tachyzoites at MOI = 3, and super- 
natants (200 )al) were collected 24 h after infection and stored at — 20°C if 
necessary. IL-12p40 levels were determined using a commercially available 
ELISA kit (Mouse IL-12 p40 ELISA kit; Thermo Fisher Scientific) accord- 
ing to the manufacturer's instructions. 

Chemokine array. For qualitative determination of the kinds of chemo- 
kines secreted by infected BMDM, chemokine Ab arrays were performed 
with culture supernatant. Macrophages were left uninfected or infected 
with VrnkuSO or VrxxkuHO Agra24 strain parasites at an MOI of 3. Culture 
supernatant was harvested after 24 h and used in mouse chemokine Ab ar- 
rays according to the manufacturer's instructions (R&D Systems). In brief, 
supernatant was incubated with an array membrane on which 25 different 
mouse chemokine Abs were bound. The membrane was then incubated 
with a mixture of biotinylated anti-chemokine Abs followed by HRP- 
conjugated streptavidin. Detection was by chemiluminescence reagents pro- 
vided by the manufacturer. 

Quantitative real-time RT-PCR. Total RNA was isolated from HFFs or 
BMDM infected with Pru or RH strains using TRIzol reagent (Invitrogen). 
cDNA was synthesized with random hexamers by using the High Capacity 
RNA-to-cDNA kit (Applied Biosystems). Samples were analyzed by real- 
time quantitative PCR for IL- 12h, IFN-y, CD70, and CD200 using TaqMan 
Gene Expression Master Mix (AppHed Biosystems) according to the manu- 
facturer's instTucdom. 2- mkroglolndin was used as an internal control gene. 

Structural studies. Based on sequence similarities, we have generated a 
GRA24-KIM 16-aa (GLLERRGVSELPPLYI) model using the HePTP 



phosphatase KIM domain (RLQERRGSNVALMLDV) from structure 
2GPH (Zhou et al., 2006) as a template, with appropriate residue mutations. 
The conformation of the side-chains and orientation of the main-chain for 
N- and C-terminal regions of peptides were modeled based on the bound 
peptide conformations in crystal structures of ERK2-pepHePTPm (2GPH) 
and p38-pepMEF2 (ILEW) peptide complexes, respectively. Similar work 
has been done with GRA24-KIM2 (5 19-534), i.e., with an R instead of 
Y, with no significant differences. A p38— GRA24-KIM binary complex was 
then prepared using Coot modeling software (Emsley et al., 2010). The 
structure of p38, isolated from the p38— pepMEF2 complex (PDB accession 
no. ILEW; Chang et al., 2002) was used as template for p38-GRA24-KIM 
docking and energy minimization calculations. The p38— GRA24-KIM 
complex was then subjected to energy minimization cycles under the Mac- 
roModel software (Schrodinger, LLC) using the OPLS-2005 (Optimized 
Potential for Liquid Simulations) force field and water as a solvent with con- 
stant dielectric value of 1. All the atoms in protein as well as peptide models 
were allowed to move freely. Energy minimization of structures were com- 
pleted via the PRCG (Polak-Ribiere Conjugate Gradient) energy mini- 
mization procedure at 10,000 iterations and performed to a convergence 
threshold of 0.02 kj mol-^-A). 

Microarray hybridization and data analysis. BMDMs were plated at 
8-10 X 10^ cells per 55 cm^ plate, infected (MOI = 3) for 18 h with KHkuSO, 
KHkuSO Agm24, FrukuSO, and VvukuSO Agm24. RNA was isolated using 
TRIzol reagent (Invitrogen), followed by phenol/chloroform/isoamyl ex- 
traction. Transcripts were obtained from three biological repHcates each. 
RNA Labehng and Array Hybridization of BMDM/RHfewSO and BMDM/ 
PrufewSO experiments were done by ArrayStar as follows: sample labehng and 
array hybridization were performed according to the One-Color Microar- 
ray-Based Gene Expression Analysis protocol (Agilent Technology). In brief, 
total RNA from each sample was linearly amplified and labeled with Cy3- 
UTRThe labeled cRNAs were purified by the RNeasy Mini kit (QIAGEN). 
The concentration and specific activity of the labeled cRNAs (pmol Cy3/)j.g 
cRNA) were measured by NanoDrop ND-1000. 1 |^g of each labeled cRNA 
was fragmented by adding 11 )jJ of lOX blocking agent and 2.2 [xl of 25X 
fragmentation buffer, the mixture was heated at 60°C for 30 min, and finally, 
55 fil of 2x GE Hybridization buffer was added to dilute the labeled cRNA. 
100 |j.l of hybridization solution was dispensed into the gasket shde and as- 
sembled to the gene expression microarray shde. The Mouse 4X 44K Gene 
Expression Array was manufactured by Agilent Technologies. The sHdes were 
incubated for 17 h at 65°C in a hybridization oven (Agilent Technologies). 
The hybridized arrays were washed, frxed, and scanned using the DNA Mi- 
croarray Scanner (part number G2505B; Agilent Technologies). Microarray 
data analysis of acquired array images was performed using the Feature Ex- 
traction software (version 10. 7. 3. l;Agilent Technologies). Quantile normal- 
ization and subsequent data processing were performed using the GeneSpring 
GX software package (version 11.5.1; Agilent Technologies). After quantile 
normalization of the raw data, genes in which at least six out of nine samples 
have flags in Detected or Not Detected (All Targets Value) were chosen for 
further data analysis. Differentially expressed genes with statistical signifi- 
cance were identified through Volcano Plot filtering between two groups 
(WT versus GRA24 mutant) of compared arrays. The thresholds are fold 
change > 2.0, p-value < 0.05. GO analysis and KEGG (Kyoto Encyclopedia 
of Genes and Genomes) pathway analysis were performed using the DAVID 
Bioinformatics Resources. Microarray data has been uploaded to GEO 
Datasets under accession no. GSE38782. 

Online supplemental material. Fig. SI shows a schematic representation 
of the strategy and analysis of GRA24 knockout. Table SI provides a list of 
the vectors and the oligonucleotides used in this study Online supple- 
mental material is available at http://www.jem.org/cgi/content/full/jem 
.20130103/DCl. 
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